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ABSTRACT 

We have carried out a large set of iV-body simulations studying the effect of residual- 
gas expulsion on the survival rate and final properties of star clusters. We have varied 
the star formation efficiency, gas expulsion timescale and strength of the external tidal 
field, obtaining a three-dimensional grid of models which can be used to predict the 
evolution of individual star clusters or whole star cluster systems by interpolating 
between our runs. The complete data of these simulations is made available on the 
Internet. 

Our simulations show that cluster sizes, bound mass fraction and velocity profile 
are strongly influenced by the details of the gas expulsion. Although star clusters can 
survive star formation efficiencies as low as 10% if the tidal field is weak and the gas 
is removed only slowly, our simulations indicate that most star clusters are destroyed 
or suffer dramatic loss of stars during the gas removal phase. Surviving clusters have 
typically expanded by a factor 3 or 4 due to gas removal, implying that star clusters 
formed more concentrated than as we see them today. Maximum expansion factors 
seen in our runs are around 10. If gas is removed on timescales smaller than the initial 
crossing time, star clusters acquire strongly radially anisotropic velocity dispersions 
outside their half-mass radii. Observed velocity profiles of star clusters can therefore 
be used as a constraint on the physics of cluster formation. 

Key words: stellar dynamics, methods: N-body simulations, galaxies: star clusters, 
star formation, open clusters and associations 



1 INTRODUCTION 

It is well known that most, if not all, stars form in star 
clusters. Star clusters can therefore be viewed as the funda- 
mental building blo cks of galaxies since the galactic disc ve- 
locity distribut ion (|Kro upa 2002, 200|}_^the galactic stellar 
mass function (jWeidner fc Kroupa 2005) and the galactic- 
field binary-star population (Kroupa 1995) are all estab- 
lished by them. 

Star clusters form as so called embedded clusters within 
the dense cores of giant molecular clouds. The star formation 
efliciency (SFE), i.e. the fraction of gas that is converted into 
stars, can be defined as follows: 



Me, 



Mecl + Mgas 



(1) 



where Med is the total mass of stars formed in the embed- 
ded cluster and Mgas the mass of the gas not converted 
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into stars. Inside molecular cloud cores, the star forma- 
tion efliciency is us ually smaller than e < 30% l|Lada|[l99^ : 
iLada fc Ladal bOOSl ). Over the whole molecular cloud, the 
star formation efli ciency is even lower and only of the order 
of a few percent (|Clark fc BonnelH 120041 ') . Since the time- 
scale for the formation of a star cluster is ~ 10® yrs, far 
larger than the typical crossing time-scale of a star clus- 
ter, star clusters are most likely in virial equilibrium be- 
fore the gas is lost, and the ve locity of the s tars reflects the 
overall gravitational potential (|Kroupal l2005) . Together with 
the low SFEs, this implies that once the primordial gas is 
expelled by UV radiation and massive stellar winds from 
OB stars or supernova explosions, star clusters will become 
super-virial and their further dynamical evolution will be 
strongly afl'ected by the gas loss. 

The impact of gas expulsion and the subsequent 
response of the star cluster has been th e subject of a 
large numb e r of theoretical investigation s jTutukovl 19781: 
Lada. Margulis fc Dearborn! I 4l I win| 

1999; Adam3 
erti 



iHilldj Il980j : 

1997al lbl: [Kroupa. Petr fc Mc Caughrean 

(2000r iKroupa. Aarseth fc Hurley! 1200 ll : ICever fc Burker 
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2OOII: iBoilv fc Kroupal l2003al lbl: iFellhauer fc Kroupal l2005l : 
Bastian fc GoodwirJ |2006| ). Using the virial theorem, 



argued that if the gas is lost instanta- 



neously and the SFE is below 50%, the entire cluster 
should be disrupted. Later analyti c and A'^-body modeling 
llLada. Margulis fc DearbornI Il984l : iGever fc BurkertI I2OOII : 
iBoilv fc KroupalbOOSal lbT have mainly confirmed Hills ideas 
but also shown that the actual limit for bound star cluster 
forma tion in case of instantaneous gas -loss is somewhat 
lower. iLada. Margulis fc DearbornI (Il984h obtained a bound 
star cluster w ith a SFE as low as e = 40%, and the 
simulations of iGever fc BurkertI (|200lh indicate that the 
critical SFE is below 35%. Using m ore realistic assumptions, 
iKroupa. Aarseth fc HurlevI (|200ll ) computed the evolution 
of Orion-Nebula-like clusters showing that they disperse 
leaving young Pleiades pre-cursors embedded in expanding 
OB associations. 

In the early phases of the evolution, before the first 
supernova explosions go off, gas is lost from star clusters 
mainly through radiation pressure and stellar winds from 
bright OB stars. Both radiation pressure and stellar winds 
need some time to expel the gas from a star cluster. A single 
10 Mq (85 M0) star for example ejects an energy of 3 • 
10^'^ (3 ■ 10^^ ) erg/s into the interstell ar medium in the form 
of radiation and mechanical energy (|Kroupall2005h . This is 
enough to unbind gas from a 1O*M0, = 1 pc cluster 
within « lO'' yrs, comparable to the crossing time of such a 
cluster. 

If the gas is not lost instantaneously but adiabatically, 
i.e. over a timescale much longer than the cluster's crossing 
time, the cluster stars can adjust to the change in potential 
and the critical SFE needed for star cluster survi v al can 
be significantly lower. iLada. Margulis fc DearbornI (119841 ) 
obtained bound clusters with SFEs as low as 20% if the 
gas expulsion timescale was equal to several initial cross- 
ing times. They also found that between 10% to 80% of 
the initial stars were lost from such systems, especially stars 
from the outer cluster parts where the initial velocities are 
more likely to exceed the escape velocities of the cluster once 
the gas is expelled. Similarly, [Goodwin, ( .1997a ) found that 
star clusters in tidal fields can form bound systems with 
SFEs as low as 20% which will survive for a Hubble time if 
the initial concentration of the cluster is high enough. Star 
clusters might also survive low SFEs if they form through 
merging of individual star clusters in star cluster complexes 
and the gas is lost in diff erent clus ters at differe nt times 
(|Fellhauer fc Kroupal lioosl ). Finally, lAdamsl (|2000l ) studied 
cluster response to gas blow out assuming that the SFE 
varies with radius, finding core survival for high central val- 
ues of the SFE. 

While many aspects of the residual gas expulsion have 
already been studied in the literature, nobody has so far 
tried to perform a systematic study of the survival limit and 
the final parameters of star clusters evolving unde r residual 
:as exp ulsio n. In addition, with the exception of iGoodwinI 
1997al ) and iKroupa. Aarseth fc HurlevI (|200ll '). the influ- 
ence of the external tidal field has so far been largely ne- 
glected. 

In the present paper, we have performed a large pa- 
rameter study of residual gas expulsion from star clusters, 
varying the star formation efficiency, the ratio of the gas 
expulsion time scale to the crossing time of the star cluster 



and the ratio of the half-mass radius to the tidal radius of 
the star cluster. Our grid of models will be useful in later 
studies of individual star clusters and also whole star cluster 
systems since the evolution of the clusters can be determined 
by interpolation between our grid points, without the need 
for further simulations. This makes it possible to determine 
the effect of gas expulsion on whole cluster systems where 
the large number of clusters prevents simulations for all in- 
dividual clusters. 

Our simulations rely on a number of simplifying as- 
sumptions, like e.g. that all stars in a cluster form at the 
same time, that the SFE is constant with radius, or that 
the overall gas density is spherical. However, relaxing most 
of these assumptions should make relatively little change to 
our results. For example, whether stars in a cluster form 
at the same time or not has no influence on the dynamical 
reaction of the cluster to gas expulsion since this reaction 
depends only on the fraction of gas thrown out. What hap- 
pens to the gas which remains, whether it has been turned 
into stars 10 Myrs before gas expulsion, or 1 Myr before, 
or at the time the other gas leaves is not relevant. Simi- 
larly, gas might leave a cluster in certain parts, where e.g. 
the overall density is lower, earlier than in other parts, lead- 
ing to a non-spherical distribution of the gas. However the 
dynamical influence of the gas is not through its density 
distribution but through the potential this density gives rise 
to and a potential field is always much smoother than the 
underlying density. Hence, the assumption that the external 
gas distribution is spherical should also have only a small 
influence on our results. 

Our paper is organised as follows: In Sec. [2] we describe 
the initial set-up of our models and in Sec.|3]we present the 
main results. SecU finally presents our conclusions. 



2 THE MODELS 

In our runs, we assumed that the SFE does not depend 
on the position inside the cluster, so gas and stars followed 
the same density distribution initially, which was given by a 
Plummer model. The gas was not simulated directly, instead 
its influence on the stars was modeled as a modiflcation to 
the equation of motion of stars. M odeling the gas as an ad - 
ditive potential has been shown bv lGever fc BurkertI (|200ll ) 
to be a good approximation of the essential physics driving 
early clust er evolution. W e used the coUisional A'^-body code 
NB0DY4 (|Aarsethlll999l ') to perform the simulations. Since 
NB0DY4 uses a Hermite scheme to integrate the motion of 
stars, the correction terms due to the external gas on the 
acceleration and its first derivative have to be evaluated at 
each step. For a Plummer model, the necessary correction 
terms can be derived analytically and are given by: 



dp 



G Mgasit) 
G Mgas{t) 



+3 



G Mgasit) 



G Mgas(t ) 

(fr) f , 



(2) 



where r is the position vector of a star relative to the clus- 
ter centre, G the gravitational constant, r = \r\, Mgas(t) 
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is the total mass in gas left at time t and rpi is the scale 
length of the Plummer model. Gas expulsion was assumed 
to start at a certai n time tp, w hich was set equal to one 
A'^-body time unit (iHeggie fc Mathicu IQSGn. equivalent to 
1/78 of a cro ssing time at the clusters virial (= gravita- 
tional) radius (|Binnev fc Tremainelll987l ). This small offset 
was introduced in order to test for each cluster if our set-up 
program created clusters that start in an equilibrium state 
when the external gas is present. After the delay time tn, 
the gas density was decreased exponentially on a character- 
istic time tm, so the total gas left at later times is given by 
l|Kroupa. Aarseth fc HurlevlbOOll ): 



(3) 



Tm will henceforth be called the gas expulsion timescale. 
The influence of the external tidal field was modeled in the 
so-called 'Near Field Approximation', which assumes that 
the size of the star cluster is much smaller than its distance 
from the galactic centre. For clusters moving in a circular 
orbit through a spherical galactic potential, the equation of 
motion of stars in a reference frame rotating with the star 
cluster can be expressed as ( Aarseth 1985,): 



dt2 



2Co X ^ -f uP'(.ixex 



(4) 



Here the first term on the right hand side is the gravitational 
acceleration due to the star cluster (both stars and gas), the 
second term is the Coriolis acceleration and the third term 
is a combination of centrifugal and tidal forces and it was 
assumed that the cluster moves in the x-y plane. The angular 
velocity to is given by: 



G Mg 



(5) 



For real star clusters, six parameters will determine 
their fate under the influence of gas expulsion: the total ini- 
tial cluster mass M^^i + Mgas, the half-mass radius, r^, and 
galactocentric distance, Rq, of the cluster, the total mass. 
Ma, of the parent galaxy inside the cluster position, and 
the star formation efficiency, e, and gas expulsion timescale, 
Tm- The number of parameters can however be reduced ef- 
ficiently: The timescale of gas expulsion for example enters 
not through its absolute value, but only through its ratio 
with the initial crossing time tcross of the star cluster. At 
the virial radius, r^, the initial crossing time tcrosa is given 
by (|Hut fc Heggiell2003l '): 



t'Cross — 2.82 



(6) 



/G ^Mgas + Med 

If gas expulsion happens slowly, stars can adjust to the 
change in potential and clusters expand adiabatically. In the 
other extreme, gas is lost instantaneously and clusters are 
strongly effected by gas expulsion. Similarly, Re and Mg 
determine together with the cluster mass the tidal radius of 
the star cluster: 

n-(^Y'Ra. (7) 



3Mg 

Dynamically important is not the absolute value of rt, but 
only the ratio of Vh/rt: If this ratio is low, clusters are nearly 
isolated and can expand freely, while if the ratio is high, they 
are strongly tidally limited and easily destroyed. 



Hence, within the framework of our model, the fate of a 
star cluster can be deduced by specifying only three param- 
eters: The star formation efficiency e, the ratio of the gas 
expulsion time scale to the crossing time of the star cluster, 
and the strength of the external tidal field, quantified by the 
ratio of the half-mass radius to the tidal radius. 

This reduction in the number of parameters makes 
it feasible to run a grid of models covering the complete 
parameter space. To this end, we have performed a set 
of A^-body simulations, varying the initial star formation 
efficiency, the ratio of half-mass radius to the tidal radius 
rh/rt, and the ratio of the gas expulsion time scale to the 
crossing time of the cluster tm /tcross- The following values 
were chosen as grid-points: 

e: 5%, 10%, 15%, 20%, 25%, 33%, 40%, 50%, 75% 

rh/n: 0.01, 0.033, 0.06, 0.1, 0.15, 0.2 
TAj/tcross-- 0.0, 0.05, 0.10, 0.33, 1.00, 3.0, 10.0 

The chosen range of parameter values essentially covers 
the relevant part of parameter space since star formation ef- 
ficiencies of less than 10% lead to cluster destruction while 
efficiencies higher than 75% have not been seen in nature 
and lead to nearly complete cluster survival. Also, clusters 
with rh/rt = 0.01 are nearly isolated and we do not expect 
that results will change much when choosing even smaller 
ratios, while clusters with rh/rt — 0.2 are strongly limited 
by the tidal field and nearly all of them are destroyed. Simi- 
larly, our values of tm /tcross cover the interesting range for 
open/globular clusters. 

All runs were perfor med with the coUisional A''-body 
cod e NB0DY4 (.AarsethI Il999l ') on the GRAPE6 comput- 
ers (iMakino et ah 20031 ) of our group at Bonn University. 
All simulated clusters contained 20.000 equal-mass stars ini- 
tially, distributed according to a Plummer sphere. The sim- 
ulations were run for 1000 initial A^'-body times (equivalent 
to about 300 initial crossing times). We found that this was 
sufficient since the gas is lost on a much shorter timescale 
and all clusters have settled into an equilibrium state by 
this time. Also, running for longer times would have meant 
that relaxation effects could have become important since 
the relaxation time of the stellar component is about 250 
initial crossing times. Our models can be viewed as a first 
step towards realistic models since two-body relaxation is 
not important in our runs and stellar evolution is neglected 
completely. In realistic star clusters, both effects would be- 
come important in the later stages of cluster evolution. 

Table 1 gives an overview of the simulations performed. 
The first three columns give the star formation efficiency 
e, the ratio of rh/rt, and the value of tm /tcross- The next 
columns give the fraction of bound mass remaining at the 
end of the runs fst = M,j /Med , the final half-mass radius in 
terms of the initial one and the global anisotropy parameter 
/?„ , which is defined by the following formula: 



Pv = 1 



(8) 



Here the sums run over all stars bound to the star cluster 
at the end of the runs and Vt and Vr are the tangen- 
tial and radial velocity component of each star. The 
anisotropies were calculated after the stellar velocities 
were transformed back into a non-rotating coordinate 
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Table 1. Details of the TV-body 
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system. For isotropic velocity dispersions, /3„ — 0, while 
radial (tangential) anisotropic velocity dispersions cor- 
respond to positive (negative) values of A complete 
list of results, including the time evolution of different 
parameters and data for runs which do not lead to cluster 
survival can be found under the following internet address: 



3 RESULTS 

3.1 Final cluster properties in dependence of 
initial parameters 



Figs. [T] and [5] depict results for the bound mass fraction 

TTT — 77 T ="^1: J — ; , ■ 1 / 7j i — J T—&t the end of the simulations as a function of star for- 

http://www.astro.uni-bonn.de/~webaiub/german/downloads.php/ 

mation efnciency, strength of the external tidal field and 

gas expulsion timescale. The different panels in Fig. [1] de- 
pict four cases for each of which the ratio of gas expulsion 
timescale to the crossing time of the clusters was held con- 



The influence of gas expulsion on the evolution of star clusters 5 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

star formation efficiency Star formation efficiency 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



Star formation efficiency Star formation efficiency 



Figure 1. Surviving mass fraction as a function of the star formation efficiency. Curves are labeled by the ratio of r^/rt used in the 
runs. The upper left panel shows runs in which the gas was removed instantaneously {tm /tcroas = 0.0). In such a case SFEs larger 
than 33% are necessary to produce bound clusters. The other panels show runs in which the gas was removed on longer timescales: 
'''M I'tcross = 0.33 (upper right), tm jicroaa = 1-00 (lower left) and Tf,.i /tcross = 3.00 (lower right). In these cases a larger number of 
stars remain bound to the clusters. In all cases, clusters in a stronger tidal field are more easily destroyed. 



stant. In case of instantaneous gas-loss (upper left panel), 
models with star formation efficiencies of 40% or larger re- 
sult in surviving clusters. In case of a very weak external 
tidal field {ru/rt = 0.01), even tlie model witfi a SFE of 
33% produces a bound cluster containing 0.9% of tlie initial 
stars, wh ich agre es will witfi estim ates from tfie literature 
ijLada. Margulis fc DearbornllTgii ). 

Prolonging the timescale on which the gas is removed 
also increases the mass fraction of bound stars in the clus- 
ters. If the tidal field is weak, star clusters can survive star 
formation efficiencies below 20% if the gas is removed slowly 



(see Table 1 and uppermost curves in the lower right panel) . 
iKroupa. Aarseth fc HurlevI (1200 ll ) performed fully realistic 
A'^-body calculations which included a tidal field realistic for 
the solar neighbourhood and gas expulsion on a thermal 
time scale. For their model B, which had TM/tcrosa = 0.3 
and rh/rt = 0.01 and a SFE of 33%, they found a bound 
mass fraction of 25% after gas expulsion, which agrees very 
well with the results of our runs. For slow gas removal, nearly 
100% of stars remain bound if the star formation efficiency 
is larger than 30%. 

Our simulations also show that the limit for survival 



6 H. Baumgardt, P. Kroupa 





Figure 2. Same as Fig. [T] but now in dependence of the gas expulsion timescale. Curves are labeled by the assumed star formation 
efficiency. Gas expulsion timescales below Tjij/tcross = 0.1 give nearly identical results. In clusters formed in strong tidal fields {ri^/rt = 
0.20), star formation efficiencies > 50% are needed in order to produce a bound system. 



depends on the strength of the external tidal field: The 
higher the ratio of Vh/rt, the higher is the star formation 
efficiency needed to produce a surviving cluster. As long 
as the strength of the external tidal field is smaller than 
rh/rt ^ 0.06, it has only a moderate influence on cluster 
dissolution. However, cluster survivability is drastically re- 
duced if the strength of the external tidal field is changed 
from rh/rt — 0.1 to Vh/rt — 0.2. In the case of rh/rt — 0.2, 
SFEs of less than 50% do not lead to the formation a bound 
star cluster since expanding clusters get easily disrupted due 
to the strong external tidal field; in this case a small ini- 
tial mass loss leads to an expansion which triggers catas- 



trophic mass-loss across the tidal boundary. Star clusters 
must therefore form concentrated in order to survive gas ex- 
pulsion. This is possible if they form for example as either 
high concentration (c > 1.5) King models or underfilling 
their tidal radius. 

Similar results can also be found in Fig. (2] which depicts 
the bound mass fraction as a function of the gas expulsion 
timescale. For gas expulsion timescales up to and including 
tm /tcross ~ 0.10, our results are nearly independent of the 
actual value of tm /tcrosa, showing that as long as the gas 
is removed with tm <^ tcross, the results do not depend on 
the details of the gas removal. Similarly, differences between 
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Figure 3. Comparison of the surviving mass fraction derived in 
this work with results from the literature. For instantaneous gas 
removal (right line of points) there is very good agreement be- 
tween both. For slow gas removal (left points, stars and open 
triangles), the critical SFE needed to produce a bound cluster 
determined here is abo ut 5% smaller than the one found by 
iGever fc BurkertI HooJ. This can be e xplained by the d i fferent 
initial density profiles and the fact that iGever fc BurkertI ll200lh 
assumed linear gas removal while we assume an exponential one. 



the TM/tcross ~ 3.0 and the tm /Tcroas ~ 10.0 case are also 
relatively small, at least for high SFEs. 

For a typical star formation efficiency of 25%, the gas 
expulsion timescale has to be larger than TM/tcross ~ 0.33 
in order to produce a bound star cluster. In our runs, star 
formation efficiencies of 5% did not lead to the formation of 
bound star clusters, which justifies that we did not perform 
any runs with lower SFE. 

Figs. [3] compares our results for the bound mass frac- 
tion as a function of the star formation efficiency with pub- 
lished results from the literature. Shown are cases when the 
gas is removed instantaneously (right group of points) and 
cases of slow gas removal. It can be seen that for instan- 
taneous gas removal there is very good agreement between 
the results of this paper and published results. SFEs of 33% 
already lead to a final bound cluster, although only a very 
small mass fraction remains bound in this case. For a SFE of 
50%, about 70% of the total cluster mass remains bound. In 
case of near adiabatic gas removal, we find that the critical 
SFE needed to produce a bound cluster is between 5% to 
10%. This is about 5% smaller than what Goyer & Burkert 
l|200lh found for their model N2 with t^xp = 10. Performing 
additional A'^-body runs shows that the difference becomes 
significantly smaller if we let the gas fraction decrease lin - 
early with time, as was done by iGever fc BurkertI (|200ll ). 
The remaining difference is probably d ue to the different 
density profiles. iGever fc BurkertI l|200ll) used King Wo = 3 
and Wo = 5 models in their runs, which are significantly less 
concentrated than the Plummer models we use. 

Fig.|4]depicts the ratio of the final half-mass radius com- 



pared to the initial one as a function of the star formation 
efficiency. In contrast to the previous two figures, only mod- 
els leading to bound clusters are plotted, which explains why 
the curves don't extend to low star formation efficiencies. 

Due to the lowering of the cluster potential when the gas 
is removed, all star clusters expand and the final half-mass 
radii are mostly larger than the initial ones. This is despite 
the fact that the tidal field removes weakly bound stars from 
the clusters. If gas is removed adiabatically, the total energy 
is conserved and the ratio of final to initial h alf-mass radius 
must obey the following relation (|Hillslll98Gl ): 



rhi 



(9) 



The dashed line in the lower right panel shows the expected 
expansion due to eq. |9| It indeed gives a good fit to our re- 
sults for nearly isolated clusters in which the gas is removed 
slowly. 

Clusters in strong tidal fields (r^/rt ^ 0.15) show only 
very small expansion as a result of the gas removal. For typi- 
cal values of the star formation efficiency («25%), expansion 
factors are around 3 or 4 in most cases. Hence, typical half- 
mass radii of embedded star clusters should be smaller by 
this factor than the corresponding ones of open or globular 
clusters. Observational data indeed suggests expansion, al- 
thou gh the ex pansion factors may be even larger than found 
here ( Kroupall2005l . The reason could be that shortly after 



gas expulsion, unbound stars are still sufficiently close to 
the cluster to be counted as members, making the clusters 
appear larger. In addition, young cluster systems contain 
a mix of bound clusters and more extended clusters which 
are dissolving but are still sufficiently compact to appear as 
star clusters. Both effects might increase the observed aver- 
age cluster radius for systems which are a few 10s o f Myr old 
(see discussion of Fig. [6l and .Bastian fc Goo dwiij (|2006l )). 

Fig. |4] also shows that in each model family, the expan- 
sion is strongest for clusters which have the smallest star 
formation efficiency that still produces a bound star cluster. 
This is to be expected since low star formation efficiency 
means a stronger decrease of the overall cluster potential 
and as long as the tidal field is weak, a cluster can expand 
since stars in the halo are not removed. For specific param- 
eter combinations, the final half-mass radius can be a factor 
10 higher than the initial one. It can also be seen from Fig. 
|4]that quick mass-loss generally leads to stronger expansion 
than adiabatic mass-loss. 

Fig. [5] finally shows the global anisotropics of the clus- 
ters by the time the simulations were stopped. Only clusters 
with e ^ 0.50 are shown since models with higher SFEs are 
hardly affected by gas expulsion. The anisotropics were cal- 
culated using eq. |5] and summing over all stars which were 
still bound to the clusters at the end of the simulations. Clus- 
ters in which the gas is removed quickly generally acquire 
radially anisotropic velocity distributions. This is due to the 
fact that such clusters are super-virial after gas removal. As 
a result, the stars expand more or less radially outwards 
and the cluster halos become populated by stars on radial 
orbits. A close inspection of these clusters shows that the 
innermost parts normally stay isotropic. Adiabatic gas-loss 
on the other hand leads to more or less isotropic velocity dis- 
persions (filled circles in Fig. [5} , at least if the clusters start 
with an isotropic velocity dispersion as is the case in our 
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Figure 4. Same as Fig.[T]for the ratio of the final half-mass radius to the initial one. The curves are labeled by the ratio of r^/rt used. 
Clusters with lower star formation efficiencies generally expand stronger than those with higher ones. Largest expansion factors are of 
order 10 and are reached for clusters which are nearly dissolved. Typical expansion factors are around 3 or 4 for SFEs of 25%, which 
means that open/globular clusters should have formed more concentrated by such a factor than as we see them now. Expansion factors 
are smaller for clusters in stronger tidal fields due to the efficient removal of stars at larger radii, and for adiabatic mass-loss. The open 
circle in the upper right panel indicates the cluster discussed in Figs.[6]to|9] The dashed line in the lower right panel shows the theoretical 
prediction for the expansion factor in case of adiabatic gas loss (eq. |9)l. 



runs. If clusters are immersed in strong tidal fields (upper 
points in Fig. [5]), the global velocity profile becomes tan- 
gentially anisotropic since the Coriolis force due to the tidal 
field forces the expanding cluster stars onto more tangential 
orbits. The clusters also acquire significant rotation in their 
outer parts in this case as a result of the non-radial accel- 
eration of the stars due to the Coriolis force, and possibly 
also angular momentum dependent escape of stars. 

The features just described might be detectable in star 



clusters which are dynamically young, i.e. for which the ratio 
of their relaxation time to their age is large tn^i/tAge ^ 1- 
A number of globular clusters in the Milky Way, like for ex- 
ample uj Cen, have relaxation times of the order of a Hubble 
time or only slightly smaller. For a subset of them, the tidal 
radius is also significantly larger than the cluster's half-mass 
radius, meaning that tidal effects are likely not important for 
these clusters. Observation of the velocity dispersion profile 
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Figure 5. Global anisotropy I3v for the final clusters with e ^ 0.50 
as a function of the gas expulsion timescale and the ratio of the 
final half-mass radius to the tidal radius. Clusters in which the 
gas was removed quickly have radially anisotropic velocity dis- 
tributions due to the sudden decrease in cluster potential (open 
circles). Adiabatic gas- loss leads to more or less isotropic veloc- 
ity dispersions (filled circles), at least if the clusters also start 
isotropic. The exception are a few clusters in strong tidal fields 
which have acquired tangentially isotropic velocity dispersions 
(triangles) due to the tidal field. These clusters are also rotat- 
ing due to the external tidal field. 



Figure 7. Evolution of projected core (lower two curves) and 
half-mass radius (upper two curves) of all stars of the two clusters 
depicted in Fig. |6] Although only the run with tm /tcross = 
0.33 results in a bound cluster, core and half-mass radius evolve 
in a very similar way in both clusters during the first 10 to 20 
crossing times, corresponding to about 10 Myr for typical cluster 
parameters (i.e. initial velocity dispersion of 5 pc/Myr). For young 
star clusters, it can therefore be difficult to distinguish bound 
from unbound clusters on the basis of their size or concentration 
alone. 



would allow to constrain the timescale over which the gas 
was expelled in these clusters. 

3.2 Evolution of individual clusters and 

observational consequences of gas expulsion 

Fig.[6]depicts the evolution of Lagrangian radii of all stars in- 
cluding unbound ones with time for two typical star clusters 
from our runs. As a result of rapid gas expulsion, both star 
clusters expand strongly in the beginning, pushing a large 
fraction of their stars over the tidal radius. For the clus- 
ter in the left panel, a small fraction of stars falls back after 
about 10 to 20 initial crossing times and forms a bound clus- 
ter. In this cluster, the Lagrangian radii become essentially 
constant with time after about 50 initial crossing times, jus- 
tifying our maximum simulation time of 350 initial crossing 
times. Compared to the initial cluster, the half-mass radius 
of the final cluster has increased by a factor 2.8 and only 
32% of the stars of the initial cluster are still bound to the 
final cluster. 

The cluster in the right panel has the same star forma- 
tion efficiency and r-^/rt ratio as the cluster in the left panel, 
but loses the gas on a slightly faster timescale {tm /tcross = 
0.10 compared to TM/tcrosa = 0.33). This small difference 
is enough to completely unbind this cluster. The dissolution 
happens within the first 50 initial crossing times. Never- 
theless, the appearance of both clusters up to this point is 
very similar in terms of the structural parameters seen by 



an observer. This is depicted in more detail in Fig. [7] which 
shows the evolution of the core radius (taken to be the 5% 
projected Lagrangian radius) and the projected half-mass 
radius of the two clusters of Fig. |6l Here again all stars in- 
cluding unbound ones were used to calculate the radii. It can 
be seen that the half-mass radii evolve in a very similar way 
initially, since even after 20 initial crossing times the half- 
mass radius of the dissolving cluster is only 40% larger than 
that of the surviving cluster. The difference in the core radii 
is larger, but is still within a factor of 4 after 20 initial cross- 
ing times. Since for typical cluster parameters, one crossing 
time is of the order of a few 0.1 Myrs, it will be difficult 
to discriminate bound from unbound clusters in young star 
cluster systems. Note also that the expansion factor of the 
projected radii is ab out 10 in bot h cases, being consistent 
with empirical data (iKroupalbOoH ). 

Fig. |8] depicts the total mass that an observer would 
derive based on the velocity dispersion, compared to the 
total mass remaining bound to the cluster. Virial masses 
were calculated for all stars projected inside the tidal radius 
of the cluster and for all stars projected inside the half-mass 
radius. The virial ratio is 3 initially since the cluster formed 
with a SFE of 33%. As the gas is removed, many stars are 
pushed over the tidal radius and the bound mass drops. 
This explains the sharp rise in the ratio of virial mass to 
the total bound mass during the first 10 crossing times. The 
virial ratio drops later as the unbound stars leave the vicinity 
of the cluster. The initial increase is weaker if only stars 
inside the half-mass radius are used to determine the virial 
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Figure 6. Evolution of Lagrangian radii of all stars, including unbound ones, in terms of the initial half-mass radius for two star clusters 
as a function of time, expressed in terms of the initial crossing time. The parameters of the clusters are given at the top of the panels. 
Solid lines show evolution of Lagrangian radii for all stars including unbound ones, dashed lines the evolution of the tidal radius of all 
stars still bound to the clusters. The cluster in the left panel expands strongly during the first 20 crossing times in response to gas loss, 
but a fraction of its stars falls back later and forms a bound star cluster. Gas expulsion happens slightly faster for the cluster in the right 
panel, which is enough to unbind it completely. 



mass due to the higher central density of the cluster and the 
smaller relative contribution of escaping stars. Nevertheless, 
even inside the half-mass radius, it takes about 30 initial 
crossing times until the virial mass estimate is a reliable 
estimate of the true cluster mass. This overestimate of a 
cluster's mass due to the early violent relaxation associated 
with gas expulsion has re cently also been pointed out by 
iBastian fc GoodwinI (|2006l ). 

Fig. |9] finally shows the anisotropy profile of the star 
cluster with e=0.33, r^/rt = 0.033 and tm /tcross = 0.33. 
The anisotropy was calculated according to eq. [S] Our ini- 
tial clusters are isotropic throughout. Due to the gas re- 
moval and the resulting escape of stars, the velocity dis- 
persion becomes strongly radially anisotropic in the outer 
parts. After 10 initial crossing times, the anisotropy reaches 
nearly /3i, — 1.0 close to the tidal radius due to escaping 
stars. After the unbound stars have left the cluster, the ra- 
dial anisotropy decreases and for the cluster shown the final 
profile is isotropic again in the center and close to the tidal 
radius, and mildly radial anisotropic at intermediate radii. 



4 CONCLUSIONS 

We have performed a large grid of simulations studying the 
impact of initial gas expulsion on the survival rate and final 
properties of star clusters, varying the star formation effi- 
ciency, ratio of gas expulsion timescale to the crossing time 
of the cluster and the strength of the external tidal field. 

Our simulations show that both the star formation effi- 
ciency and the speed with which the gas is removed have a 



strong influence on the evolution of star clusters. In the case 
of instantaneous gas removal, clusters have to form with 
SFEs ^ 33% in order to survive gas expulsion. This limit 
is significantly lowered for gas removal on longer timescales 
and clusters with SFEs as low as 10% can survive gas expul- 
sion in the adiabatic limit if the external tidal field is weak. 
External tidal fields have a significant influence on the clus- 
ter evolution only if the ratio of ru/rt is larger than about 
0.05. Below this value, star clusters behave nearly as if they 
are isolated. 

All star clusters expand due to gas expulsion. For star 
formation efficiencies around 25%, star clusters expand by a 
factor of 3 or 4, nearly independent of the speed with which 
the gas is removed. Open or globular clusters must therefore 
have formed more compact and with higher central densities 
than with what we see them today. Star clusters that just 
managed to hold together should have expanded strongest 
and in rare cases expansion factors of 10 or more are possible. 

The velocity dispersion profiles of surviving clusters are 
tangentially anisotropic in their outer parts if the clusters 
are strongly effected by the external tidal field. For nearly 
isolated clusters, the velocity profile is radially anisotropic if 
the gas was removed on a timescale smaller than the crossing 
time of the cluster, while slow gas expulsion leads to mainly 
isotropic velocity dispersions. In the central parts of star 
clusters, the velocity profile preserves its initial shape. This 
offers a chance to test how star formation and gas expulsion 
is happening in nature by studying the velocity profiles of 
sufficiently isolated and dynamically young {tnei > tAge) 
star clusters. 
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Figure 8. Mass derived from the observed velocity dispersion 
compared to the true mass remaining bound as a function of 
time for the surviving cluster of Figs. |6] and [7] Shown are mass 
estimates based on all stars inside the tidal radius (solid line) 
and inside the cluster's half-mass radius (dashed line). Due to 
the assumed SFE of 33%, the initial total mass (stars plus gas) 
is three times higher than the stellar mass alone. The estimated 
mass increases compared to the total mass due to gas removal 
and the loss of a large fraction of the cluster mass. Depending 
on which stars are used to derive the total cluster mass, it takes 
between 30 to 50 initial crossing times (i.e. about 20 Myrs) until 
the virial mass becomes a reliable estimate for the total mass. 

The simulations reported here should be useful for a 
number of follow-up projects. First, they allow us to study 
the impact of gas expulsion on the evolution of the mass 
function and other properties of whole globular cluster sys- 
tems in galaxies: Starting with a range of cluster parameters 
and SFEs, the mass fraction remaining bound to each indi- 
vidual cluster can be calculated by interpolating between 
the runs in our grid. Repeating the process for all clusters 
will then give the impact of gas expulsion on the whole star 
cluster system and as a function of the age of the system. 
Our runs also allow to test to which degree mass determi- 
nations of young star clusters through measurements of the 
radial v elocity dispersion are affec ted by the initial gas ex- 
pulsion l|Bastian fc Goodwinll200^ . Finally, observations of 
individual star clusters can be directly compared with our 
simulations to infer their starting conditions. We plan to 
carry out such projects in the near future. 
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Figure 9. Anisotropy parameter /3„ as a function of radius (ex- 
pressed by the enclosed mass fraction) at three different times for 
the same cluster as in Fig. [8] The initial cluster model (solid line) 
is isotropic throughout. In the beginning, the cluster becomes 
strongly radial in the outer parts due to gas removal which puts 
many stars on outward-bound, ballistic orbits (dotted line). The 
final profile at the end of the run is nearly isotropic in the center 
and close to the tidal radius and mildly radially anisotropic at 
intermediate radii (dashed lines). 
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